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Abstract 

We report the first observation of tlie flavor-changing neutral current decay B K*i'^i^ and 
an improved measurement of the decay B Ki^i~ , where i represents an electron or a muon, 
with a data sample of 140 fb~^ accumulated at the T{AS) resonance with the Belle detector at 
KEKB. The results for the branching fractions are B{B K*e+i-) = (11.5l|2±0.8±0.2) x 10"'^ 
and B{B Ki+i-) = {4:.8tli ± 0.3 ± 0.1) x 10"^, where the first error is statistical, the second 
is systematic and the third is from model dependence. 

PACS numbers: 11.30.Hv, 13.20.He, 14.65.Fy, 14.40.Nd 
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Flavor- changing neutral current (FCNC) processes are forbidden at tree level in the Stan- 
dard Model (SM); they only proceed at a low rate via higher-order loop diagrams. SM decay 
amplitudes for the FCNC processes B — >■ Xs^ and B Xs£^i~ , where Xs denotes inclusive 
hadronic final states with a strangeness S = ±1 and i represents an electron or a muon, 
have been calculated with small errors . If additional diagrams with non-SM particles con- 
tribute to these FCNC processes, their amplitudes will interfere with the SM amplitudes, 
making these processes ideal places to search for new physics 

Measurements of the decay rate for B — > Xg'j as well as the recent first exclusive 
and inclusive measurements by Belle for B Ki'^i~ j3| and B Xsi'^i~ [5] have so far 
shown no disagreement with the SM predictions. Deviations due to non-SM amplitudes 
are often expressed in terms of the Wilson coefficients C7, Cg and Ciq; a strong constraint 
on the magnitude of C7 has been set by i? — >• X^'j, and a large area of the C9-C10 plane 
has been excluded by i? ^ Ki~^£~ and B — > Xsi~^£~ 6]. A complete determination of all 
three Wilson coefficients, including the sign of C7, requires the measurement of the forward- 
backward asymmetry in B ^ K*i~^i~ or B —>■ Xsi~^i~] however, B — > K*i~^i~ has not been 
previously observed. 

In this Letter, we report the first observation of the decay B K*i^i~, using a data 
sample of 152 million B meson pairs, corresponding to 140 fb~^ taken at the T{AS) reso- 
nance. We also report an improved measurement oi B ^ K£^£~ , superseding our previous 
result based on 29 fb"^ 0. 

The data are produced in e^e~ annihilation at the KEKB energy-asymmetric (3.5 on 8 
GeV) collider and collected with the Belle detector j^. The Belle detector is a large- 
solid-angle spectrometer that includes a three-layer silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel threshold Cherenkov counters (ACC), time- 
of-flight (TOF) scintillation counters, and an electromagnetic calorimeter (ECL) comprised 
of CsI(Tl) crystals located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron fiux-return located outside of the coil is instrumented to identify 
muons (KLM). 

The event reconstruction procedure is similar to our previous report j^. We reconstruct 
the following final states: 5° ^ K*^i+i-, B+ K*+i+i-, 5° ^ iTg^+f" and B+ 
K^i^i~ . Charge conjugate modes are implied throughout this Letter. The following decay 
chains are used to reconstruct the intermediate states: K*^ K'^tc', K*^ —>■ K^tc'^ and 
K*+ K+Ti^, K^g TT+TT-, and 7r° ^ 77. 

Charged tracks are classified as e, /i, K and vr candidates by discriminating between the 
fiavors for the pairwise combinations, using criteria which allow multiple classifications of an 
individual track. The e/h discriminant (where /i = or vr) is formed from the energy deposit 
in the ECL, the specific ionization measurements in the CDC, and the ACC light yield. The 
discriminant is based on the hits in the KLM. The K/n and K/^ discriminants use 
the CDC, ACC, and TOF information. Electrons, muons and kaons are selected using loose 
conditions on the e/h, fi/h and K/tt discriminants, respectively. All tracks are classified 
as pions unless they satisfy tight conditions on e/h or K/tt; the same e/h condition is 
required for kaons, and a similar K/fi condition is required for muons. To reduce the 
misidentification of hadrons as leptons, we require minimum momenta of 0.4 GeV/c and 
0.7 GeV/c for electrons and muons, respectively. We apply a tight requirement for the 
muons below 1.0 GeV/c. Each of the charged tracks, except for the Kg tt+tt^ daughters, 
is required to have an impact parameter with respect to the interaction point of less than 
0.5 cm transverse to, and 5.0 cm along the positron beam axis. Photons are reconstructed 
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within the ECL with a minimum energy requirement of 50 MeV. 

Invariant masses for the vr", and K* candidates are required to be within windows 
of ±10 MeV/c2 (~2cr), ±15 MeV/c^ (~3.3a) and ±75 MeV/c^ (1.51), respectively, around 
their nominal masses. We require a minimum momentum of 0.1 GeV/c for the vr'' candidates. 
We impose selection criteria based on the distance and the direction of the K'g vertex 
and the impact parameters of daughter tracks. For K*'^ K~^7i^, cos6'hei < 0.8 is required 
to reduce background from soft vr^s, where ^hei is the angle between the K*~^ momentum in 
the B rest frame and the momentum in the K*~^ rest frame. 

We form B candidates by combining a K^*^ candidate and an oppositely charged lepton 

pair using two variables: the beam-energy constrained mass Mbc = \J {E^ca.m/ ~ \v*b/^\'^ 
and the energy difference AE" = — E^^^^, where p*^ and E'^ are the measured momentum 
and energy, respectively, of the B candidate, and -E^eam is the beam energy. Here and 
throughout this Letter, variables denoted with an asterisk are calculated in the T{4S) rest 
frame. When multiple candidates are found in an event, we select the candidate with the 
smallest value of |Ai?|. 

The following five types of backgrounds are considered. 1) Charmonium B decay back- 
ground from B J/ipXs and B —>■ ip'Xg decays is removed by vetoing lepton pairs whose 
invariant mass (M«) is near the J/ip or ip' mass Q. In addition, we reject events that have 
a photon with energy less than 500 MeV within a 50 mrad cone around either the electron 
or positron direction (or a photon within each cone) and an e+e~7(7) invariant mass within 
the veto windows. For K*i~^i~ modes, we reject the event if an unobserved photon along 
one of the lepton directions with an energy -Ebeam ~ ~ replace the pion, giving 

Mu^ and Mbc consistent with B ^ J/ip K. 2) We suppress background from photon conver- 
sions and 77° e"'"e~7 by requiring the dielectron mass to satisfy M^e > 0.14 GeV/c^. This 
eliminates possible background from B K*'y and B K^*\^ . 3) Background from con- 
tinuum qq (g = u, d, s, c) production is suppressed using a likelihood ratio T^cont formed from 
a Fisher discriminant, cos^^, and, for K^*^e^e~ only, cos^gp^. The Fisher discriminant 
is calculated from the energy flow in 9 cones along the B candidate sphericity axis and the 
normalized second Fox- Wolfram moment R2 The angles 9^ and 9*^^ are the B meson 
angles with respect to the beam and the sphericity axes, respectively. 4) Semileptonic B 
decay background is suppressed using another likelihood ratio T^sb formed from the missing 
energy of the event, -Ej^iss; cos9*^. 5) Hadronic B decay background, B K^*^h^h~, 
e.g., from B Dtt, can contribute if two hadrons are misidentified as leptons. We find that 
other potential backgrounds are negligible. 

For each decay mode, the selection criteria on the two likelihood ratios T^cont and T^si are 
chosen to maximize A^5 / A/iV5~±~/VB , where Ns is the expected signal yield and Nb is the 
expected background in the Mbc and AE signal windows. The signal windows (~2.5cr) are 
defined as |Mbc - Mb\ < 0.007 GeV/c^ for both lepton modes and -0.055(-0.035) GeV < 
AE < 0.035 GeV for the electron (muon) mode. A large Monte Carlo (MC) background 
sample of a mixture of 6 — c decays and e"'"e"" — > qq events is used to estimate Nb- The 
B — s> K^*H~^i~ signal events are generated according to Ref. to determine Ns, and to 
estimate the efficiencies that are summarized in Table IH 

The signal yield is determined by a binned maximum-likelihood fit to the Mbc distribution 
for the events within the AE signal window using a Gaussian signal plus three background 
functions. The area of this Gaussian function is the signal yield; the mean and width are 
determined using observed J/ ip K^*^ events. We find no dilepton mass dependence of the 
width and mean using a MC study. The first background function is for the semileptonic B 
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decays and, to a lesser extent, the continuum background, and is modeled with a threshold 
function 111] whose shape parameter is determined using a large MC sample that contains 
oppositely charged leptons and whose normalization is allowed to float. This MC sample 
reproduces the background parametrization for B K^*^e^iJi^ data in which only combina- 
torial background is expected. The two other background functions account for the residual 
B to charmonium decays and hadronic B decays, and are modeled with separate combi- 
nations of a similar threshold function and an additional Gaussian component. The shape 
and the size of the charmonium background function are fixed from J/ if) and if)' inclusive 
MC samples. We find the Gaussian component of the charmonium background contributes 
less than one event. The shape and the size of the hadronic background are evaluated 
using hadron enriched data by relaxing the lepton identification criteria. The Gaussian 
components of the hadronic background contribution, multiplied by the lepton misidenti- 
fication probability (measured in bins of momentum and polar angle with respect to the 
positron beam), are then found to be 1.05 ± 0.08 and 0.64 ± 0.05 events for B Ki^i~ 
and B — >■ K*i~^i~ , respectively. 

Figure [T] and Table IJ give the fit results. We observe 35.8t7;3(stat.)±1.7(syst.) B —>■ 
K*£~^£~ signal events with a significance of 5.7, and 37.9l6;9(stat.)^J'°(syst.) B K£^£~ 
signal events with a significance of 7.4. The error due to uncertainty in the fixed parame- 
ters is included in the systematic error. To evaluate the uncertainty in the signal function 
parametrization, the mean and width of the Gaussian function are changed by ±1 stan- 
dard deviation (a) from the values determined from J/ipK^*^ events. The uncertainty in 
the semileptonic plus continuum background parametrization, which is the largest error 
source, is obtained by varying the parameter by ±lcr from the value determined with a large 
MC sample. The uncertainties of the hadronic (charmonium) background contributions are 
evaluated by changing the shape parameters and the normalizations of the Gaussian and 
threshold components by ilex (±100%). The significance is defined as y'— 2 ln(£o/'Cmax), 
where £max is the maximum likelihood in the Mbc fit and Co is the likelihood of the best 
fit when the signal yield is constrained to be zero. In order to include the effect of sys- 
tematic error in the significance calculation, we use the parameters simultaneously changed 
by la (100% for the charmonium background) in the direction that reduces the resulting 
significance. 

In addition to the systematic error in the signal yield, we consider the following ex- 
perimental systematic errors in the efficiency determination. For each charged track, we 
estimate the systematic error due to reconstruction efficiency to be 1.0%, and the system- 
atic errors due to kaon, pion, electron and muon identification to be 1.0%, 0.8%, 0.5% and 
1.2%, respectively. For each candidate and candidate, we estimate the systematic 
errors due to reconstruction efficiencies to be 4.5% and 2.7%, respectively. The uncertainty 
in the background suppression is estimated to be 2.3% using J/ijjK^*^ control samples. Sys- 
tematic errors due to MC statistics range from 0.5% to 2.2%. All these errors are added in 
quadrature. 

The uncertainty due to the theoretical model assumptions is evaluated by calculating the 
efficiency for signal MC samples generated using three form-factor models jy, 0| and taking 
the maximum difference as the model-dependence error. 

When calculating the branching fractions, we assume an equal production rate for charged 
and neutral B meson pairs, isospin invariance, lepton universality for B — > K£~^£~, and the 
branching ratio B{B K*e^e~)/B{B K*^^^~) = 1.33 6]. The combined efficiency and 
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branching fraction are scaled to the muon mode. We find 

B{B K*tr) = (11.5^1^ ± 0.8 ± 0.2) x 10"^ 
B{B Ktr ) = (4.8+J:°'± 0.3 ± 0.1) X 10"^ 

where the first error is statistical, the second is systematic, and the third is from model 
dependence. This systematic error is a quadratic sum of the systematic errors in the yield 
and efficiency, and the uncertainty in B meson pair counting of 0.5 %. The results are 
consistent with the SM predictions [1, [13, , our previous values , and results recently 
reported by BaBar . The complete set of results is given in Table E] 

For the modes with a significance of less than 3, we set 90% confidence level upper limits. 
The upper limit on the yield, A^, is defined as /J^ C{n)dn = 0.9 /q°° C{n)dn. The function 
C{n) is the likelihood for signal yield n, using signal and background shape parameters that 
are modified by la of their errors in the direction to increase the signal yield. The upper 
limits for the branching fractions are then calculated by using the efficiencies reduced by la 
of their errors. 

Figure 121 shows the measured = Mf^c^ distributions for B Ki'^i" and K*£^i~. The 
signal yield is extracted in each bin from a fit to the Mbc distributions. 

In summary, we have observed the decay B K*i~^i~ for the first time. This mode will 
provide a useful sample for a forward-backward asymmetry measurement. The B — > Ki^i~ 
decay is also measured with improved accuracy. The measured branching fractions are in 
agreement with the SM predictions, and may be used to provide more stringent constraints 
on physics beyond the SM. 
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FIG. 1: Mbc distributions (histograms) for K^*'>l'^l samples. Solid and dotted curves show the 
results of the fits and the background contributions, respectively. 
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FIG. 2: distributions of (a) K*i~^£~' and (b) Ki^i^ . Points with error bars show the data while 
the hatched boxes show the range of SM expectations from various models 0, 1^ . Statistical and 
systematic errors are added in quadrature. 
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